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The ratio of the B + B~ and B°B° production rates in e + e~ annihilation is computed as a function 
of the B meson velocity and BB*n coupling constant, using a non-relativistic effective field theory. 
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The dominant production mechanism for B mesons at 
CLEO, BaBar and Belle is via the P-wave decay of the 
T(4S) state, e+er -> T(45) -> BB. The final state 
can contain either charged (B + B~) or neutral (B°B Q ) 
mesons, and the ratio of charged to neutral B mesons 
produced enters many B decay analyses, including stud- 
ies of CP violation. We will define the ratio 



R + '° = 1 + 5R + '° 



r (T(4g) -> B+B~) 
T (T(45) B°B°) ' 



(1) 



which is unity in the absence of isospin violation. The ex- 
perimental value measured by the BaBar Collaboration 
is R+/° = 1.1O±O.O6±O.O50, and by the CLEO collabo- 
ration is 1.04±0.07±0.04 gfand 1.058±0.084±0.136 @. 

Isospin violation is due to electromagnetic interactions, 
and due to the mass difference of the u and d quarks. In 
most cases, isospin violation is at the level of a few per- 
cent. However, it is possible that there can be significant 
isospin violation in T decay J§ §, §. The T(45) is barely 
above BB threshold; the B mesons are produced with a 
momentum ps ~ 338 MeV and velocity v/c — 0.064 [us- 
ing M T(4S) = 10.58 GeV, M B = 5.2792 GeV], so that the 
final state is non-relativistic. The electromagnetic contri- 
bution to is a function of v and the fine-structure 
constant a. In the non-relativistic limit, there are l/v 
enhancements, and the leading contribution is a function 
of a/v, 



i?+/° = 



ira/v 



or 
Av 2 



ira _ / a 
= 1+ 2^+°'- 



,(2) 



and can be obtained by solving the Schrodinger equa- 
tion in a Coulomb potential for a P-wave final state . 
Corrections to this result are suppressed by powers of 
a without any l/v enhancements. For T(AS) decay, 
this gives P + /° = 1.19,a significant enhancement of the 
charged/neutral ratio ||[ [| ||. 

Lepage [pi computed corrections to Eq. (^) by assum- 
ing a form-factor at the meson-photon vertex, and found 
that <5P + /° could be significantly reduced from 0.19, or 
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FIG. 1: One-loop correction to T(4S I ) 
ton and (b) pion exchange. 
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even change sign. Recent advances in the study of heavy 
quark systems and non-relativistic bound states allow us 
to improve on this estimate of <5P + /°. Since the final 
state B mesons are non-relativistic, and have low mo- 
mentum, the final state interactions of the B meson can 
be treated using non-relativistic field theory combined 
with chiral perturbation theory || . At momentum trans- 
fers smaller than the scale of chiral symmetry breaking 
A x ~ 1 GeV ||, the photon vertex can be treated as 
pointlike. The B and B* states have a mass splitting 
of 45.78 ± 0.35 MeV, which is small compared with the 
momentum of the B meson, so the B and B* must 
both be included in the effective theory. Since is 
much smaller than the mass of the 6-quark, heavy quark 
spin symmetry holds |l0|, pjj , and one can treat the B 
and B* as one multiplet described by the field of 
HQET |1. Similarly, the B and B* can be combined 
into a field, whose properties are related to those 
of by charge conjugation |ll|. At low velocities, 
the dominant isospin violation is that enhanced by fac- 
tors of l/v, which is obtained by solving the Schrodinger 
equation with the — interaction potential. The 
NRQCD counting rules |Q show that BB annnihilation 
is suppressed, and can be neglected. At low momentum 
transfer, the H^ b ' — H^ b ' potential is dominated by single- 
pion exchange. Isospin violation in the potential arises 
from Coulomb photon exchange, and from isospin viola- 
tion in the pion sector due to the 7r + — tt° mass difference 
and r] — ir° mixing. 

In perturbation theory, the first contribution to 5R + /° 
is from the graphs in Fig. [j]. The one-loop photon graph 
gives the ira/2v term in Eq. (||). It is enhanced by tt 2 /v 
compared with a typical relativistic radiative correction, 
which is of order cx/tt, because of the non-relativistic na- 
ture of the integral. The one-loop pion graph is similarly 
enhanced by ir 2 /v ~ 150 compared with a typical chiral 
loop correction. As a result, the correction from Fig. 
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is not small, and cannot be treated in perturbation the- 
ory. However, it is possible to sum the multiple pion 
exchanges by solving the Schrodinger equation using the 
one-pion plus one-photon exchange potential. This sums 
the series of graphs shown in Fig. |^. Additional correc- 
tions, such as vertex corrections, are not included in the 
Schrodinger equation. However, these corrections are not 
enhanced by it 2 /v, and so are subleading compared with 
the terms we have retained. 

The H^—H^' interaction potential is the same as the 
jjw — fjv>) potential (by charge conjugation symmetry), 
and was computed in Ref. |l5| which studied bbqq exotic 
states. The potential depends on the B*Bir coupling 
constant g which is not known. Heavy quark symmetry 
implies that g is the same as the D*Dtt coupling. The 
D* can decay into Dir (via the coupling g) or D7 (via 
electromagnetic interactions), and the decay rates can 
be used to obtain g A fit to the experimental 

data gives two possible solutions, g = 0.27 ^02 -o'o2 or 
g = 0.76 to.ol -OA |H, with the smaller value being pre- 
ferred. A recent measurement of the D* + widt h by the 
CLEO collaboration gives g = 0.59 ± 0.01 ± 0.07 |9j. We 
will give our results as a function of g. 

The T(4S I ) is a 1" state, and can decay into five 
possible channels, (i) BB with S = 0, I = 1, (ii) B*B* 
with S = 0, i = 1, (iii) B*B* with S = 2, I = 1, (iv) 
B*B* with S = 2, 1 = 3 and (v) BB* +B*B with 5=1, 
1=1, where t is the orbital angular momentum and S is 
the total spin. Since the T(45) is below BB* and B*B* 
threshold, only the first state is allowed as a final state, 
but all five states need to be included as intermediate 
states in the calculation. [The actual number of states 
is double this, since one has both charged and neutral 
channels.] Let 77, /3 = 1 — 5 denote one of the five possible 
£S states and a,b — 1,2 denote the charged and neutral 
sectors, respectively, so that a given channel is labeled by 
the index pairs r\a or (5b. The radial Schrodinger equation 
has the potential 



where £ v — (1,1,1,3,1) are the angular momenta of the 
various channels. The denominator is uib since is 
the reduced mass of the BB. The Coulomb potential is 



V 1 1 



r) 



(6) 



where a is the fine-structure constant. It only contributes 
to the charged sector a = b = 1, and does not mix differ- 
ent tS states. 

The pion potential can be computed using the tech- 
niques given in Ref. ||, EL ||: 



Kja,/3b( r ) — 



2U rl p(m„+,r) hlU n p (%o,r) 



a b 



where h\ = 1.01, = 0.99 [|2j. The structure of the 
potential is easy to understand. The off-diagonal ele- 
ments are transition amplitudes between the charged and 
neutral sectors due to 7r + exchange, and depend on the 
charged pion mass vn^ and coupling constant g, which 
is included in the definition of U. The diagonal matrix 
elements are due to tt° exchange, and depend on m n a. In 
the absence of r\ — ir° mixing, the ir + coupling constant 
is y/2 time the n° coupling which gives the 2 : 1 ratio of 
the off-diagonal to diagonal elements. The values of h+ 
and ho differ from unity due to 77 — ir° mixing |2^] . 

The computation of the matrix U ri p{m. n ,r) is non- 
trivial. The answer is that 



where 



U(m, r) 



U(m, r) = TU (m, r) T 



g 2 m 2 e mr 
8irf 2 r 



(8) 
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V]a,d r ) + KM r ) + Mv^a.b (3) 



where V* is the pion potential, V 1 is the Coulomb po- 
tential, V is the angular momentum potential, and M v 
is the contribution due to the B* — B mass difference, 
Am. 



ui(m, r) — H 

mr 



u 2 (m,r) 



f ~ 132 MeV is the pion decay constant and 



(9) 



Mi = 0, M 5 = Am, 
M 2 = M s = Mi = 2Am. 



(4) 



The B° - B+ mass difference is 0.33 ± 0.28 MeV @| , and 
will be neglected in our analysis. Note that a B° — B + 
mass difference of 0.33 MeV contributes about 0.05 to 
5R+/° from the p 3 dependence of the phase space of the 
P-wave decay. 

The angular momentum potential is 



Vr>a,8b( r ) 



'/ 



1) 



-S a ,bSr)fd 



(5) 
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and T* is the transpose of T. 

Equation (|^) is the leading contribution to the long 
distance part of the potential. As argued in Ref. JL5], 



FIG. 2: Series of graphs summed by solving the Schrodinger equation. The dashed line represents pions and photons. 



Eq. (Q) will dominate the potential until r ~ 1/(2771^) at 
which point two-pion exchange begins to contribute. We 
introduce a cutoff r m ; n = l/(2m,r), and use Eq. (Q) for 
r > r min , and set V""' = for r < r min . The short distance 
part of the potential can be included into a renormaliza- 
tion of the production vertex. The Coulomb potential 
will be allowed to act until r = 0. 

The T(45) is produced by the space component of the 
electromagnetic current tr/ l b. Heavy quark spin symme- 
try holds in the T system |l(| juj, so the T(45) decays 
into H0>) - such that the spins of the heavy quarks 
in the final mesons are combined to form the spin of the 
T(4iS), i.e. the polarization of the virtual photon. The 
orbital angular momentum and spin of the light degrees 
of freedom are combined to form total angular momen- 
tum zero. A little Clebsch-Gordan algebra shows that the 
amplitude for the T(AS) to decay into the five channels 
is 

The amplitude for decay to the I — 3 channel is zero 
to this order in the velocity expansion. The coefficients 
c a , a = 1,2 are unknown, but the absolute values of c a 
are irrelevant for the computation of i? + /°; all that is 
needed is the ratio c\jcb of the charged to neutral pro- 
duction amplitudes. The dominant production of the B 
mesons is via the isosinglet T(45) state, in which case 
ci = C2- Isospin violating effects, including direct pro- 
duction of B's not via the T(45) lead to a deviation of 
Cx/c% from unity. As discussed above, cutoff effects in the 
potential can be absorbed into the production amplitudes 
c a . One expects short-distance corrections to introduce 
isospin violation in the ratio c\jci of a few percent, the 
typical size of other isospin violating effects in hadron 
physics. We will define 5c by c\ /c2 = 1 + Sc. The value of 
5c is related to the value of r min , since changes in the cut- 
off induce changes in the Lagrangian coefficients. Since 
5c is unknown, our computation of is uncertain at 

the 5% level; however the uncertainity is much smaller 
than the expectation that 5R + /° is 19% from Coulomb 
interactions alone. Cutting off the Coulomb potential at 
short distances reduces the value of 5R + ^°. Since the 
Coulomb potential is the dominant source of isospin vio- 
lation, one expects that 5c will be negative. 

The method of computation is as follows. One solves 
the Schrodinger equation with potential Eq. (0). The 
boundary condition on the wavefunction as r — > oo is 
that one has a plane wave plus an outgoing scattered 
wave. [One can see this directly from the sum of graphs in 
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FIG. 3: R +/0 as a function of g for v = 0.064 and 5c = 0.02 
(dotted), 0.0 (solid), -0.02 (dashed) and -0.04 (dot-dashed). 



Fig. |.] Only the B+B- and B°B° states exist as prop- 
agating modes as r — > oo; the other channels have expo- 
nentially decaying wavefunctions. The plane wave state 
is chosen to be in the B + B~ or B Q B° channels to com- 
pute the charged or neutral meson production rates, re- 
spectively. The overlap of the computed wavefunction as 
r — ► with the production amplitude Eq. (O) gives the 
final production amplitude, the absolute square of which 
gives the production rate. [Note that the wavefunction 
near r = can have all five channels.] The answer for 
depends on 5c, g and the velocity v of the outgoing 
B meson. Provided the dominant production mechanism 
is via the photon coupling to the heavy quark, the result 
for R+/° holds even away from the T(4S) resonance since 
it depends only on the quarks being non-relativistic. The 
value of c a will depend strongly on the beam energy, and 
peak at the resonance, but 5c, the isospin violation in 
the production amplitude should be a smooth function 
of energy. 

In Fig. U we have plotted as a function of g for 

5c and v = 0.064, the value in T(45) decay. R +/0 is 
approximately constant and equal to its value from only 
Coulomb corrections, Eq. (0), until g > 0.6, at which 
point starts to decrease. is approximately 

constant for small g even though the shifts in the pro- 
duction amplitudes are large. The one-loop pion graph 
in Fig. [j] is about three times the tree-level graph. Sum- 
ming the pion graphs in Fig. || gives about a 20% (for 
g ~ 0.6) shift in the charged and neutral production 
rates. The rates into the charged and neutral channels 
vary by about a factor of two for the range of Yukawa 
couplings in Fig. |l but their ratio varies by about 

10%. For larger values of g than those shown, has 
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FIG. 4: R +/0 as a function of v for g = 0.3 and 5c = 0.02 
(dotted), 0.0 (solid), -0.02 (dashed) and -0.04 (dot-dashed). 
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FIG. 5: R +/0 as a function of v for g = 0.8 and 5c = 0.02 
(dotted), 0.0 (solid), -0.02 (dashed) and -0.04 (dot-dashed). 



rapid v dependence due to the formation of meson bound 
states, because the pion-exchange potential is sufficiently 
attractive. For our choice of parameters, this occurs for 
g ~ 1.3, well outside the allowed range EJ, |i"sf] . 

In Fig. U and ||, we have plotted as a function 

of velocity for different values of 5c for two illustrative 
choices g — 0.3 and g — 0.8 consistent with the two 
solutions for g found in Ref. . The vertical line is the 
velocity at the T(45). At the T(45) peak, for g = 0.8, 
varies from 1.17 to about 1.09, whereas for g = 0.3, 
varies between about 1.25 and 1.1. 

In Figs. ^| and [?|, we have plotted as a function 

of v for g = 0.3 and g = 0.8, respectively, for different 
values of the cutoff from r m j n = 1/(2771^) to For 
small values of g, the variation of the cutoff does not 
change For larger values of g, the cutoff variation 

is consistent with expectations from naive dimensional 
analysis || . A factor of two variation in the cutoff intro- 
duces a 4% variation in R + /°. 

The absolute value of R + /° depends on the value of 
5c, and the cutoff r nim . If g is small (~ 0.3, the pre- 
ferred value in Ref. jl8|), then for values of 5c consistent 
with expectations from dimensional analysis, one expects 
<5i? + /° > 0.1. The Yukawa corrections do not signih- 
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FIG. 6: R + ^° as a function of v for g = 0.3, 5c = 0, and cutoffs 
rmin = l/(2m,r) (solid) , l/V2m-K (dashed), 1/mn (dotted). 
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FIG. 7: R +/0 as a function of v for g = 0.8, 5c = 0, and cutoffs 
r m ; n = l/{2m 7r ) (solid), l/vlm^ (dashed), \jm-n (dotted). 



cantly change from the Coulomb value. We note, 

however, that this is due to a cancellation in R + /° after 
summing the graphs in Fig. ^ the one loop pion correc- 
tion from Fig. [3] is about three, and is not small. If g 
is close to the larger value g = 0.8, then 5R + ^° at the 
T(45) is smaller, but still around 0.1. In this case, there 
is some cutoff dependence, so is more uncertain. 

The dependence of on v (or equivalently, y/s) is 

calculable. One can see that the curves in Fig. Q have a 
different shape from those in Fig. [|, so measuring R + I a 
as a function of v can provide information on the B* Bit 
and D*Dtt coupling g, which is needed for many calcula- 
tions, such as the ratio of the B s — B s to B — B mixing 
amplitudes Jl3| . 
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